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Abstract 

The time-varying fluid and optical fields of several 
cavity configurations have been computed on over- 
set mesh systems using the Reynolds-averaged Navier- 
Stokes equations and geometric optics. Comparisons 
between numerical results and Airborne Optical Ad- 
junct (AOA) flight data are made in two-dimensions 
for a quieted cavity geometry with two lip-blowing 
rates. In three-dimensions, two proposed aero-window 
locations for the Stratospheric Observatory For In- 
frared Astronomy (SOFIA) are discussed. The sim- 
ulations indicate that convection of large shear layer 
structures across the aperture cause the blur circle di- 
ameter to be three times the diffraction-limited diam- 
eter in the near-infrared band. 

Nomenclature 

c speed of sound 

dB decibel, 20 log 10 

/ frequency 

h enthalpy 

I intensity 

k wave number, ^ 

K ratio of convection by freestream speed 

L characteristic length 

m mass flow rate 

M Mach number 

MTF modulation transfer function 

n index of refraction 

OPD optical path difference 

p instantaneous static pressure 

q velocity magnitude or dynamic pressure 

Re Reynolds number 

St Strouhal number, 

SR Strehl ratio, = exp( — $ 2 ) 

t time 
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T absolute temperature 

T c characteristic time, — 

u, r, w Cartesian velocity components 

x^y,z Cartesian physical space or 

aperture coordinates 
A angstrom, 10“ 10 m 

a angle of attack 

3 Gladstone-Dale constant, (n — I)srp 

0 momentum thickness 

A wavelength 

p density 

a shear layer spreading rate parameter 

phase, 2 jlOPD 

( ) mean quantity 

< > root mean square quantity 

( )' fluctuating quantity, / = / + / / 

Subscripts 

STP standard temperature and pressure 

T total quantity 

o© freestream quantity 

Introduction 

The study of light propagating through an unsteady 
fluid field has important applications ranging from 
laser weaponry to astronomy platforms. Airborne 
housing of these systems provides mobility, mainte- 
nance, and performance advantages which, in combi- 
nation, can be superior to land or space-based alter- 
natives. However, prediction of the fluid and optical 
behavior of these airborne systems remains a difficult 
problem. 

This report describes the progress of a computa- 
tional approach for use in the design of transonic 
aero-windows. The prediction methodology has been 
driven by the design of the Stratospheric Observatory 
For Infrared Astronomy (SOFIA), the successor to the 
Kuiper Airborne Observatory (KAO), which will offer 
ten times the resolution of the KAO. Figure 1 depicts 
the SOFIA, which will have a 2.5 meter Cassegrain 
telescope mounted in a cavity of a Boeing 747SP. In 
order to numerically assess the safety and performance 
of this platform, extensive evaluation of the computa- 
tional methods by comparison against experiment is 
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Fig. 1: Artist’s concept of the SOFIA configuration 


necessary. 

Many studies of the effect of a fluid field upon 
an optical field have been conducted over the past 
four decades. Many experimental and theoretical ap- 
proaches to the optical distortion problem have been 
investigated; only those of which are pertinent to this 
transonic aero-window problem are summarized here. 
The experimental efforts can be grouped into two cat- 
egories: direct measurement methods and techniques 
based on aerodynamically inferred quantities. Results 
obtained via the latter method are more prevalent be- 
cause of practicaJ difficulties in direct measurement 
techniques. 1 In fact, only aerodynamically inferred dis- 
tortion levels will used for validation of the present 
work. 

Although early experimental and theoretical efforts 
assumed incoherent statistical turbulence, 2 3 4 recent 
studies have begun to examine the effect of shear layer 
structures on electromagnetic field distortion. Using 
a passive scalar field from a direct numerical simu- 
lation, Truman and Lee 5 found an optimum viewing 
angle normal to the hairpin vortices in the homoge- 
neous sheared fluid region. They also found analysis 
via non-refracting geometric optics to be equivalent to 
the parabolized Helmholtz representation of light. Al- 
though this class of studies provides excellent insight 
into the effects of small-scale structure on the electro- 
magnetic field, it is clear that the computational ex- 
pense of such methods precludes their near-term use 
for the problems under consideration here. 

The study of large scale structures in shear layers 
has been an active topic of research since they were 
observed by Brown and Roshko in 1974. 6 Only re- 
cently has the effect of these structures on the optical 
field been studied. In 1990, Chew and Christiansen 7 8 
experimentally observed the effect of shear layer struc- 
tures on beam propagation. Tsai and Christiansen 9 
used an Euler simulation to determine the optical char- 
acteristics of a perturbed free shear layer. It was hy- 
pothesized that the effect of the vortical structures on 
the optical field could be modelled by a growing sinu- 
soidal phase plate. 

The numerical modelling of the optical effect of a 
cavity-spanning shear layer was presented by Cassady, 
Birch, and Terry 10 in 1987. They found their two- 


dimensional solution to result in poor prediction of 
optical distortion. Farris and Clark 11 12 used time- 
mean quantities and empirical evidence to ascertain 
the fluctuating density levels required for optical anal- 
ysis. 

The present effort attempts to determine what por- 
tion of the optical path distortion can be resolved using 
cell sizes required to obtain an accurate flowfield so- 
lution. Towards this end. computed optical distortion 
levels are compared to flight or wind tunnel measure- 
ments for two- and three-dimensional quieted cavities. 

Previous reports 13,14 have described the method 
development for two-dimensional free shear layers, a 
backward-facing step, and a rectangular cavity. 15 Com- 
parison of the computed cavity case with Rossiter's 
data showed agreement in the dominant resonant 
peaks to within 5 dB. The computed and experimen- 
tally observed pressure loading trends were similar 
along the cavity walls. In three-dimensions, rectangu- 
lar and treated quiet cavity solutions were computed 
and compared to experiment. 16 Sound pressure lev- 
els along the cavity bulkheads for both the resonating 
and quieted geometries were found to be in agreement. 
However, although the power spectra of the experi- 
ment and computation were similar at low frequen- 
cies, numerical dissipation caused a rapid decrease in 
energy content at high frequencies. 

Although the optical model has been described pre- 
viously, 13 this paper documents the extension of the 
model and provides new validation information. The 
following sections address the method used to predict 
the unsteady flows and the resultant optical distortion. 
Analysis of the aperture fluid and optical fields for 
AOA and SOFIA configurations are presented. For the 
two-dimensional AOA geometry, time- varying density 
fields and optical path lengths are shown. Short and 
long exposure far-field diffraction patterns are com- 
puted for a three-dimensional aft cavity SOFIA con- 
cept. 

Approach 

Solution of the aircraft and cavity fiowfields were 
computed using models for the fluid field, the effect 
of turbulence, and the optical distortion. A diagonal 
scheme 1 ' was used for the solution of the Reynolds- 
Averaged Navier-Stokes equations, implemented in an 
overset grid framework. 18 Euler implicit time integra- 
tion and second-order spatial differencing was used, 
with viscous impermeable wall conditions specified as 
no-slip, zero normal pressure gradient, and adiabatic. 
Information transfer across overset mesh boundaries 
was implemented using trilinear interpolation of the 
dependent variable vector, Q — [p s pu, pv, p\L\e] T . Al- 
gebraic turbulence models were used, implemented 
with a variable F maz cutoff for wall-bounded flows and 
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a shear layer model for the cavity aperture region . 19 20 
The flow solver cost was 13// .s/cell/iteration on a single 
head of the Numerical Aerodynamic Simulator (NAS) 
Cray Y/MP-832. 

Generally, a significant effort in grid generation 
is required before flow and optical analysis can be- 
gin. However, recent advances in algebraic 21 and 
hyperbolic 22 methods have enabled rapid discretiza- 
tion of complex geometries. Hyperbolic grid genera- 
tion, which provides spacing and orthogonality control, 
was used for the wall- bounded regions, while algebraic 
grids were used in shear flow regions including plumes 
and wakes. This choice of topology allows simple spec- 
ification of turbulent regions and also permits the recy- 
cling of meshes, useful for configuration changes such 
as cavity positioning. 

The optical computations documented here use a 
refracting-ray method, reported on earlier, which is 
limited to studying the effects of the resolved large- 
scale structures. 23 The method tesselates a structured 
grid into tetrahedra and uses piecewise mean indices 
of refraction for each of the tetrahedra. Indices of re- 
fraction were computed using n = 1 + P — , where 
the Gladstone-Dale constant, /?, can be found using 
the Cauchy formula. 

Assessment of the optical performance of an aero- 
window begins by specification of the ray initialization 
plane. Integration of the optical path length through 
the aero-window is then performed along the rays at 
specified time increments in both the streamwise, x, 
and crossflow, y, aperture directions. The resultant 
OPD{t % x,y) can be used in a complex aperture func- 
tion of the form P = A Assuming no loss in 

transmittance, then the wave amplitude A is unity in 
the aperture and zero elsewhere, while the phase of the 
wave is computed from OPD(t<x,y) = The au- 
tocorrelation of the complex aperture function, P * P, 
gives the far field diffraction pattern, computed using 
a two-dimensional Fourier transform. Time averag- 
ing successive short-period diffraction patterns gives a 
long exposure result. Integration of the intensity of 
this resultant long or short exposure diffraction pat- 
tern gives the area for a specified encircled energy level. 
From this area the equivalent blur-circle diameter due 
to the resolved fluid scales can be found. Inclusion of 
the effects of small-scale turbulence could be incorpo- 
rated into the computation of long exposure blur cir- 
cles by multiplication of the above modulation transfer 
function (MTF) with a turbulence MTF. 3, 24 

Results and Discussion 

Aero-optical simulations of the U.S. Army Airborne 
Optical Adjunct (AOA) and the SOFIA configurations 
are discussed below. Information pertaining to the 
computation and analysis of the unsteady flowfields. 


including grid resolution and turbulence modelling, is 
given elsewhere. 14, 23 


2-D AOA Cavities 

Data available from flight tests of the AOA 23 shown 
in Fig. 2, provides valuable validation information for 



Fig. 2: U.S. Army Airborne Optical Adjunct 


the present simulations. These two-dimensional nu- 
merical simulations were used to determine if optical 
quieting methods, particularly aft ramp treatment and 
lip-blowing, could be accurately simulated. The flow 
about the geometry, depicted in Fig. 3, was computed 
in conjunction with two lip-blowing rates for the for- 
ward aperture only. The 1009c lip-blowing rate case 
corresponded to a m = 0.42(/3u)oo. For the discussion 
below, computed high and low lip-blowing rates refer 
to 1009? and 1%; of this mass flow rate. 

M = ,75 


.1 

0 




Fig. 3: AOA case: instantaneous Mach contours 

Computed and flight mean Mach number profiles 
are compared in Fig. 4 for two lip-blowing rates, with 
instantaneous Mach number contours shown above 
each set of profiles. The quantity d> indicates the angle 
from the cupola crest at which the data was measured. 
Overall, the Mach number contours, which were aver- 
aged over 2000 time steps, show agreement in the max- 
imum vorticity as a function of x-station. However, 
differences exist between experimental and computa- 
tional results at the low-speed edge of the shear layer. 
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Fig. 4: AOA case: instantaneous Mach number contours 
and mean profiles at (a) low and (b) high lip-blowing rate 




Fig. 5: AOA case: power spectra, low lip-blowing rate 


which may be due to blockage in the cavity of the air- 
craft that was not computationally represented. The 
discrepancy at the upper aft portion of the shear layer 
appears to be due to blockage in the computational 
model. 

Comparisons of pressure spectra at the aft ramp for 
the low lip-blowing rate are shown in Fig. 5. The com- 
puted spectra can be seen to be quantitatively and 
even qualitatively different from flight dp” a. The com- 
puted result lies more than 15 dB below ae data, and 
a peak in the low lip-blowing rate spectra is clearly 
computed, but is not seen in the flight data. The high 
lip-blowing rate spectra were similar to that shown in 
Fig. 5, albeit without the spectral peak at 340 Hz. 23 

It has been noted from experimental evidence 26 that 
the frequency of large structures in shear layers is inde- 
pendent of axial station and occurs at Strouhal num- 
ber of St = = 0.024 ± 0.003. where 9 is the 

local shear layer momentum thickness and / denotes 
frequency. This phenomena is corroborated bv the re- 
duction of other researchers' data, 6 7 9 ■ 2l ' 28 whose 
results range from about St = 0.02 to 0.03 for incom- 
pressible shear layers. 

Momentum thickness can be estimated using 
Gortler's solution, giving 6 = 0.036|^7-r for <to = 11.0, 
which compares favorably to the empirically deter- 
mined correlation 29 of 9 = 0.034j-^x. Using this re- 
lationship along with a compressibility correction, 30 
the computed peak in the AOA solution at 340 Hz 


corresponds to a Strouhal number of 0.032. For com- 
parison, peaks were found in SOFIA experiments and 
computations at approximately St = 0.028. 

Based upon these experimental observations, it is 
hypothesized that large scale shear layer structures are 
being resolved. However, the lack of empirical sup- 
port from the flight data pressure power spectra is at 
odds with this conclusion. The comparison is further 
clouded by the reasonable comparison in < p' > for the 
low lip-blowing rate shown in Fig. 6, and the presence 
of the organized structures shown in Fig. 7. 
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Density fluctuation, 

Fig. 6: AOA case: Sp > profiles with (a) low and (b) 
high lip-blowing 
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Three-dimensional effects are a possible explana- 
tion for the discrepancy between computation and 
flight. Rockwell 31 noted that for sufficiently large 
Reynolds numbers three-dimensionality reduces coher- 
ence in the shear layer. This implies that an error in 
the assumption of planar flow for the small flow os- 
cillations considered in these quieted cavity configu- 
rations. The evolution of streamwise-oriented vortic- 
ity interacting with the primary vortices would act to 
spread peaks in the reattachment ramp pressure spec- 
tra. As a final note. Fig. 5 also depicts data, the ordi- 
nate scaled by (?oc ) flight /(<hc) tunnel and the abcissa 
by (Coc/ L) flight/ (coc./ L)tunneh obtained from an AOA 
wind tunnel test. 32 The wind tunnel data shows that 
a small peak exists where expected according to the 
above analysis. 



Fig. 7: AOA case, low lip-blowing rate: instantaneous 
contours of (a) and (b) with schematic of optical 
model 

Using the computed unsteady density field, aero- 
optical effects can be determined. Figure 6 compares 
the computed and experimental profiles of the root 
mean square of the density fluctuations. Levels of 
< p f > were computed over a time segment of about 
90 ms in increments of 0.44 ms. Using the elapsed 
time for a particle to convect across the aperture at the 
mean shear layer speed as a characteristic time, then 
the optical computation was taken for about nine T c . 
In Fig. 6, y is the rake angle from horizontal, with the 
axis of rotation offset from the cupola centerline. De- 
termination of the systematic error band on the exper- 
imental result is discussed below. The low lip-blowing 
rate result underpredicts the magnitude of the peak 
in > , however the peak location is in fair agree- 
ment. The computed results for the high lip-blowing 
rate compare poorly to experiment, possibly due to 
inadequate grid resolution and/or the increased flow 
complexity of the merging shear layers. This type of 
active control is presently not a design option for the 
SOFIA, therefore further effort toward improvement 
of the high lip-blowing case was not warranted. 

Further investigation of the low-blowing rate case 
revealed the presence of large convecting structures as- 


sociated with t lie shear layer. Figure 7 shows a contour 
plot of , depicting the growth and propagation of 
these shear layer structures. Also depicted in Fig. 7 is 
a schematic of the optical model, with the initial and 
final stations of the optical path integration given by r 0 
and rf. The large structures, associated with a 0.03^ 
vertical velocity component, are the primary contribu- 
tors of the computed density fluctuations of the shear 
layer. The speed of the waves, as determined from 
Fig. 8, is 0. 561/30, below the value of 0.661/30 inferred 
by Rossiter 10 for rectangular cutouts, yet above the 
0.51 determined analytically by Roscoe and Han- 
key. 33 
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Fig. 8: AOA case: contours of OPD f (xJ)[ in.] along 
aperture, low lip-blowing rate 

In 1990, Chew and Christiansen 8 and Tsai and 
Christiansen 9 deduced that a free shear layer model 
of a sinusoidal phase delay growing in x would pro- 
duce results similar to those observed in both compu- 
tation and experiment. Figure 8 displays behavior of a 
similar nature for the aero- window problem of concern 
here. 

Comparisons of integrated aero-optical quantities, 
shown in Fig. 9, reveal slight overprediction for the 
low lip-blowing case and, given the < p f > profiles, 
expected underprediction for the high lip-blowing rate 
case. Also shown in Fig. 9 are the < OPD f > for two 
additional integration paths, shown to demonstrate 
that most of the optical distortion is caused by the 
shear layer. 

The result for the integration path which extends 
from r 0 = -8" to y = 12" displays an increment in 
< OPD ( > of about (7 x 10~ 7 )" from the 7" path 
length case which passes through the shear-layer alone. 
The path initialized above the shear layer, from r 0 = 
4" to Tf = 12", shows a small < OPD' > indicating 
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that < p > a bow the shear layer is small. Finally, 
the time mean optical path difference, OPD, ran f>e 
seen to contribute curvature to the wavefronts a* the 
light propagates through the shear layer. The optical 
clarity of the shear layer was determined using a J = 
2.584 x l(r 4 , matching the value which was used to 
reduce the experimental data. 



Fig. 9: AOA case: streamwise variation of optical path 
difference quantities 

The analytic result for the < OPD f >. which goes 
like j\ is found from 34 

OPD' >\ 2 



Derivation of the model, which utilizes time-mean 
quantities to determine < OPD f >, is given by 
Bogdanoff. 34 This analytic result assumes an index- 
matched shear layer with a sinusoidal n profile, n(y) = 

ni ~ - r m “ sin ( T7~ — V The constants, — = 0.0091 
and Tson.m ~ 1.316 w = 1.31(0.18x[m]). are empirical 



relations. The virtual origin of the shear layer is placed 
ar i = 0" to obtain this analytic result. As shown in 
Fig. 9, the analytic gradient in < OPD f > is in good 
agreement with flight data. 


Reduction of Experimental Data 

The reduction of the data obtained from experi- 
ment 30 is noted here to delineate the approximations 
used and estimate systematic error bounds in the op- 
tical path distortion levels. Values of p' are computed 
from assumptions of quasi-steady flow: 


h T = c p T (l + V— 
Differentiation with respect to t gives 


RT f r = 



+ uu' 


7 ~ 1 

7 


Using {pu) f = pu' 4- p'u then 


5 = - + (7 - 1)A/ 2 — - [(7 - DA/ 2 + 1] p - 

1 p pu 1 1 p 


The experimental observations against which the com- 
puted results are compared assume simultaneously 
small fluctuations in pressure and total tempera- 
ture, 36, 3 ‘ resulting in 


t 

P 


f— ^ 

\(l - DA/ 



' (Pu)' 
pu 


(1) 


Mean Mach number and density profiles are deter- 

/ w 

mined from isentropic relations, while is propor- 
tional to the voltage fluctuation, obtained from 
hot film probes. The optical path disturbance is then 
found from 1 


(OPD 1 ) 2 = 2(-P) f L < p' > 2 l r dr (2) 

V PST P ) Jo 

where l -£ is the turbulent eddy size relative to the shear 
layer width, determined from cross correlation data to 
be typically about 15%. 

The few available independent measurements 37 38 
indicate that pressure fluctuations of about 2% of 
freestream static pressure occur in the shear layer 
spanning the aperture of a quieted cavity geometry. 
In fact, Hahn 25 reported pressure fluctuations of 8% 
from shear layer rake measurements, however these in- 
clude the dynamic pressure component normal to the 
orifice as well. Pressure fluctuation levels can also be 
inferred from sound pressure levels in the cavity, ob- 
served to be at least 130 dB for the AOA case. Shear 
layer total temperature fluctuations of about 1%. have 
also been reported for this Mach regime 38 The present 
low lip-blowing computation found a < p f >% 1% and 
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a < Tj O.8/1 in the shear layer. The assumption 
of < V T >, < // >w 0 in a shear layer is therefore 
questionable, and is used to estimate systematic ex- 
perimental error bounds. 

The determination of the error in < p' > due to 
background noise levels begins by assuming the pas- 
sage of a compression wave parallel to the static pres- 
sure port in the wake rake. Normal reflection of the 
wave would impart a larger deviation f: a < p f > 

as computed by Eq. 1. Utili zing Gbrtler s free shear 
layer solution to provide u(r), assuming a cavity tem- 
perature recovery factor of unity, and holding mean 
static pressure constant through the layer, then p(r) 
is defined. The sensitivities of ^ to ^ and ^ are 

P p TV 

± (7 -iiAf ; + i and respectively. Using 

a compression or rarefaction wave of strength < p' > 
through the shear layer, then local values of p 1 due to 
wave passage are defined. This value of p f provides 
the error bound about the value obtained from Eq. 1, 
which assumes negligible < p* > and < Tj >. Taking 
shear layer pressure fluctuation levels corresponding 
to 135 dB and a velocity ratio r = 0.1, then the sys- 
tematic error in the density fluctuations is 0.13% at 
the shear layer center. Figure 6 shows the resultant 
systematic error bars in < p f >. 

From Eq. 2 the value of < p 1 > is linearly propor- 
tional to < OPD f >. The error in < OPD 1 > can be 
found by using a conservative within-system error of 
0.05% in > gleaned from Fig. 6. plus the systematic 
error from tfie above analysis. The resultant error bar 
i plotted in Fig. 9. 



Fig. 10: Forward cavity SOFIA cases: Comparison of 
sound pressure levels 

namic measurements were used to infer distortion. 

Comparison of experimental and numerical density 
fluctuations are shown in Fig. 11, where < p' > can be 
seen to be severely underpredicted. Although peaks 
in the density fluctuations were computed, the highly- 
ordered shear layer structures similar to those found in 
the AOA study were not observed. Differences may be 
attributable to grid coarseness or within-system mea- 
surements errors. 
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Forward Cavity SOFIA 

Wind tunnel tests, completed in 1990, allowed vali- 
dation of cavity acoustic response and optical charac- 
teristics. The cavity environment results are summa- 
rized in Fig. 10 for both resonant and quieted config- 
urations. 13 Aerodynamic measurements in the shear 
layer were used to infer optical quantities for the qui- 
eted cavity, configuration 100. 

Following computation of the unsteady flow, the op- 
tics code was applied to the computed density field 
obtained for configuration 100 from f = 0 to 7.8 ms. 
Ten rays were propagated through 110 instantaneous 
density fields in time intervals of At = 70.6 ps. The op- 
tical measurement was taken for about five again 
using the elapsed time for a shear layer structure to 
convect across the aperture as the characteristic time. 
The forward cavity SOFIA results presented here are 
for a computational plane at approximately the cross 
flow center of the aperture, which will provide only a 
streamwise variation in optical properties. The numer- 
ical results are presented compared to previous anal- 
ysis 34 and experiment 16 in which shear layer aerody- 






■ 


2 

1 

0 

• 1 

.... Exp. (Roc* A Cooley) 
— Present 

HH Estimated error bend 

V' 

1 X-37.8” 

\ X«41.$” 


p- 

[ x»45.r* 

< 

.01 .02 .030 .01 02 .03 0 

.01 02 03 


Denelty fluctuation. 


Fig. 11: Configuration 100: density fluctuation at cross- 
flow center of aperture 


The optical wavefront distortion through the con- 
figuration 100 aero-window is summarized in Fig. 12. 
Figure 12a shows that the distortion model applied 
through the shear layer alone underpredicts the data 
determined analytically and experimentally. However, 
the computed trend is generally consistent with the 
data. At the streamwise center of the aperture, the 
aerodynamically inferred < OPD f > at two additional 


spanwise locations arc shown. These points provide an 
estimate of the crossflow variation in experimental dis- 
tortion levels. 
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Fig. 12: Configuration 100: comparison of wavefront 
distortion 

Figure 12b depicts computed < OPD 1 > for ray 
propagation originating below the secondary mirror, 
r 0 = —3.7", and above the shear layer, r 0 = 2.3". 
Comparison of the computed results show an incre- 
ment in < OPD ' > below the secondary mirror. This 
distortion increment appears to be caused by a jet of 
re-entrant fluid originating from the shear layer im- 
pinging on the aft ramp. Finally, Fig. 12c shows that 
curvature is imparted to the mean optical field. The 
dip in the fluctuating and mean OPD levels at x = 42" 
is caused by the presence of the secondary mirror, in 
which the index of refraction, n, was fixed at unity. 

Aft Cavity SOFIA 

Forward placement of the telescope in a favorable 
pressure gradient region has an advantage in terms of 
an optically thin boundary layer. However, the fuse- 
lage moldline and structural complexities forward of 


the wing present considerable manufacturing difficul- 
ties. An alternative site for the telescope aft of the 
wing reduces the modification costs and permits the 
use of a larger usable cavity volume. However, an aft 
cavity site has potential problems of scattered light 
omitted from engines and plumes, an optically thick 
boundary-layer, unknown cavity response, and possi- 
bly poor empennage flow behavior at off-design condi- 
tions. 

Figure 13 depicts the simplified geometry used to 
address some of these concerns: horizontal and vertical 
stabilizer geometry was unavailable for this simulation. 
Details of this flight condition simulation are available 
elsewhere. 14 



Fig. 13: Aft SOFIA case: Surface C p and plume tem- 
perature contours 


The acoustic response of the aft cavity is compared 
to scaled data from the forward cavity experiment in 
Fig. 14. The computed result is taken from a location 



Fig. 14: Aft SOFIA case: Power spectra comparison 

on the aft ramp, while the experiment is from a loca- 
tion within the cavity. The agreement of spectra is rea- 
sonable to about 100 Hz, above which grid coarseness 
dissipates energy rapidly. 14 Figure 14 shows peaks at 
60 and 110 Hz, the latter corresponding to a Strouhal 
number of 0.028. 

During the aft cavity SOFIA computation the en- 
tire aperture density field was saved in increments of 
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0.G8 m.s\ every five flow solution steps. Using this den- 
sity field, propagation of a plane wave through the 
aperture revealed variations in the wavefront distor- 
tion. as shown in Fig. 13. These ordered variations 
in OPD . indicative of shear layer structures in the 
aperture, impact the aft ramp at a frequency of 110 
Hz. giving a St — 0.028. Figure 15 shows maximum 
distortions of about one wavelength. ±A q, with a re- 
sultant maximum < OPD 1 > of approximately 0.7 Xp. 
Computation of the OPD(t. j\ y) was performed for a 
64 x 64 array of rays normal to the aperture and ini- 
tialized just above the secondary mirror. The optical 
integration was performed over 8 T c . 



Fig. 15: Aft SOFIA case: Sample wavefront distortion 
history 


Using these phase distortion levels, far field diffrac- 
tion patterns were then computed. Figure 16 depicts 
the diffraction-limited Airy pattern for reference, and 
both instantaneous and time-averaged exposures. The 
instantaneous exposure pattern shows some evidence 
of speckle, with a large reduction in central peak in- 
tensity. This spreading of energy is manifested in the 
computed Strehl ratio of 0.34. 

Finally, Fig. 17 shows that the large scale structures 
in the shear layer cause the equivalent 80% blur cir- 
cle to be three times the diameter of the diffraction- 
limited case. However, as can be seen in Fig. 16, the 
blurring in the streamwise direction is worse than in 
the crossflow direction. Note that the because the 
small scale fluid motion is modelled when using the 
Reynolds-averaged Navier-Stokes equations, the com- 
puted blur circle is much smaller than actually ob- 
served. 39 


2.5 m aperture 
X = X D = 5893 A 



instantaneous exposure 



28 ms exposure SR=0.34 


0 12 3 

Arcsec 

Fig. 16: Aft SOFIA case: Far field diffraction patterns 
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Fig. 17: Aft SOFIA case: Normalized integrated intensity 
distributions 

Conclusions 

Computations of quieted cavity configurations have 
shown convection of large scale flow structures across 
the aperture. The shedding frequency of these struc- 
tures compare reasonably well with experimentally de- 
termined shear layer Strouhal numbers. The com- 
puted results indicate that three-dimensional effects on 
the shear layer spanning a quieted cavity can be signifi- 
cant. The differences in two- and three-dimensional re- 
sults are manifested in the spectra of the pressure loads 
and in the magnitude of the optical wavefront distor- 
tion. Since the primary contributors to the computed 
OPD' w'ere the large scale structures, computations of 
the Strehl ratio were found to be reasonable. However, 
because the small scale fluid motion is modelled, the 
blur circle diameter is significantly underpredicted. 

Further improvements to the prediction of optical 
performance may be found from investigation of shear 
regions with direct Navier-Stokes methods 5 or use of a 
turbulence MTF. 24 Finally, although low Reynolds and 
Mach number number experiments show large struc- 
ture formation, direct OPD measurements at realistic 
conditions would be useful for validation of numerical 
aero-optical studies of this type. 
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